Comparison of renormalization group schemes for sine-Gordon type models 
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The scheme-dependence of the renormahzation group (RG) flow has been investigated in the 
local potential approximation for two-dimensional periodic, sine-Gordon type field-theoric models 
discussing the applicability of various functional RG methods in detail. It was shown that scheme- 
independent determination of such physical parameters is possible as the critical frequency (temper- 
ature) at which Kosterlitz-Thouless-Berezinskii type phase transition takes place in the sine-Gordon 
and the layered sine-Gordon models, and the critical ratio characterizing the Ising type phase tran- 
sition of the massive sine-Gordon model. For the latter case the Maxwell construction represents 
a strong constraint on the RG flow which results in a scheme-independent infrared value for the 
critical ratio. For the massive sine-Gordon model also the shrinking of the domain of the phase 
with spontaneously broken periodicity is shown to take place due to the quantum fluctuations. 

PACS numbers: ll.lO.Gh, ll.lO.Hi, ll.lO.Kk 
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I. INTRODUCTION 

Since the invention of the renormalization group (RG) 
method it has been the main goal of functional 
RG to describe systems where the usual approxima- 
tions (e.g. perturbation theory) failed. Strongly cor- 
related electrons, critical phenomena, quark confinement 
are examples where the non-perturbative treatment is 
required. The RG equations arc functional integro- 
differential equations, consequently, they can only be 
handled by using truncations (for reviews see The 
truncated RG flow depends on the particular choice of 
the so-called regulator function of the RG method, i.e. on 
the renormalization scheme, hence, the predicting power 
of the RG method is weakened. Indeed, once approxi- 
mations are used a dependence of physical results on the 
choice of the regulator function is observed, consequently, 
any RG scheme makes sense if this dependence is weak. 
Therefore, it is of relevance to clarify how far the results 
obtained are independent (or at least weakly dependent) 
of the particular choice of the renormalization scheme 
used. 

The purpose of this paper is to investigate the RG 
scheme-dependence of two-dimensional (2D) periodic 
scalar field theories with possible inclusion of explicit 
symmetry breaking mass terms, i.e. the sine-Gordon 
(SG) the massive sine-Gordon (MSG) [i| and the 
multi-component layered sine- Gordon (LSG) Q mod- 
els. Our motivation is twofold: (i) from the methodi- 
cal point of view, 2D SG-type models represent a new 
platform to consider RG scheme-dependence since the 
comparison of results obtained by various RG schemes 
has been investigated previously for the 0(N) symmetric 
polynomial scalar field theory mostly in three-dimensions 
B S Hj S [lo| (ii) from the phenomenological point of 
view, 2D SG-type models have important direct real- 
izations in high-energy and condensed matter physics, 
consequently, their physical parameters should be deter- 



mined independently of the particular choice of the RG 
scheme used. 

In particular, our goal is to consider under which con- 
ditions it is possible to determine scheme-independently 
three physical parameters, the critical frequency of 
the SG theory, the critical ratio u/M^ of the MSG model 
and the layer-dependent critical value 01{Nl) of the LSG 
model. The SG scalar field theory (the Bose form of the 
massive Thirring fermionic model) represents the sim- 
plest non-trivial quantum field theory which is used to 
study the confinement mechanism. It has two phases 
separated by the critical value of the frequency 01 — Stt. 
The 2D SG model belongs to the universality class of 
the 2D Coulomb gas and the 2D-XY spin model which 
have been applied to describe the Kosterlitz-Thouless- 
Berezinskii (KTB) phase transition of many con- 
densed matter systems. Due to the exact mapping of the 
generating functional of the SG model onto the partition 
function of the Coulomb gas, 01 is related to the criti- 
cal temperature of the condensed matter system [l^, HI] 
which can be measured directly. Moreover, the MSG 
model describes the vortex dynamics of 2D supercon- 
ducting thin films [3, [lB| and it is the bosonized version 
of the one-flavor massive Schwinger model, i.e. the 2D 
quantum electrodynamics (QED2) .4, ^IfJ, [3 which 
possesses confinement properties. The Ising type phase 
transition of the MSG model is controlled by the di- 
mcnsionlcss quantity u/M"^ where u is the Fourier am- 
plitude and M is the scalar mass. It is related to the 
critical ratio of QED2 which separates the confining and 
the half- asymptotic phases of the fermionic model which 
has been calculated by semiclassical and lattice meth- 
ods [H, [13, The LSG model is used to describe 
the vortex dynamics of magnetically coupled layered su- 
perconductors [1^ m, m, m, HHI and considered to be 
the bosonized form of the multi-flavor Schwinger model 
(i.e. multiflavor QED2) [lEIMIlil. It is known that the 
critical frequency of the LSG model (and consequently 
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the critical temperature of the corresponding condensed 
matter system) depends on the number of layers Nj^ in a 
particular way PUNl) = %ttNlI{Nl~1) [15]. This laj^r 
number dependence can be observed experimentally 27j. 
Therefore, the critical frequency 01 of the SG theory, 
the critical ratio u/AP of the MSG model and the layer- 
dependent critical value P'^{Nl) of the LSG model are 
physical parameters, consequently, it is an important is- 
sue to consider under which conditions it is possible to 
determine them scheme independently which is the main 
goal of the present work. 

The structure of our paper is the following. Some gen- 
eral remarks on RG scheme-dependence is given in Sec- 
tion ini The applicability of various RG schemes to the 
2D SG model and the scheme-dependence of the criti- 
cal frequency is discussed in Section IIIII In Section 
IIVI we investigate the scheme-dependence of the critical 
ratio of the MSG model which characterizes the bound- 
ary of the phases of the model. In Section |V] we con- 
sider the dependence of the critical frequency (31{Nl) 
of the LSG model on the RG scheme used. Section 
IVII presents the summary and our concluding remarks. 
A brief overview of the frequently used RG methods is 
given in Appendix [X] in order to remind the reader on 
their advantages and drawbacks and settling the nota- 
tions used throughout the paper. Finally, in Appendix [Bl 
the scheme-dependence of the renormalization of the one- 
component polynomial scalar field theory is reviewed. 



II. MOTIVATION AND GENERAL REMARKS 
ON SCHEME-DEPENDENCE 

In the framework of the Kadanoff- Wilson [l[ RG ap- 
proach the differential RG transformations are realized 
via a blocking construction, the successive elimination of 
the degrees of freedom which lie above the running ul- 
traviolet (UV) momentum cutoff k. Consequently, the 
effective theory defined by the (bare or effective) blocked 
action contains quantum fluctuations whose frequencies 
are smaller than the momentum cutoff. This procedure 
generates, e.g. the functional RG flow equation. 
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for the effective action Tk[<f\ when various types of reg- 
ulator functions Rk are used, where T), ' [4>] denotes the 
second functional derivative of the effective action (see 
e.g. [2l). Here Rk is a properly chosen infrared (IR) 
regulator function which fulfills a few basic constraints 
to ensure that Vk approaches the bare action in the UV 
limit (fc — > A) and the full quantum effective action in 
the IR limit (fc 0). Indeed, various renormalization 
schemes are constructed in such a manner that the RG 
flow starts at the bare action and provides the effective 
action in the IR limit, so that the physical predictions 
(e.g. fixed points, critical exponents) are independent of 
the renormalization scheme particularly used. 



Since RG equations are functional partial differential 
equations it is not possible to solve them in general, 
hence, approximations are required. One of the com- 
monly used systematic approximation is the truncated 
derivative expansion where the (bare or effective) action 
is expanded in powers of the derivative of the field, 



Vk{(i^)- 



Zk{<i>)\{d^<i>f 



In the local potential approximation (LPA) higher deriva- 
tive terms are neglected and the wave-function renormal- 
ization is set equal to constant, i.e. Zk = 1. The solution 
of the RG equations sometimes requires further approx- 
imations, e.g. the potential can be expanded in powers 
of the field variable (with a truncation at the power 
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where the scale-dependence is encoded in the coupling 
constants Cn{k). Since the approximated RG flow de- 
pends on the choice of the regulator function the physical 
results could become scheme-dependent. 

A general issue is the comparison of results obtained by 
various RG schemes [1, H B S US] which has been investi- 
gated for the 0(N) symmetric polynomial scalar fleld the- 
ory in three dimensions in great detail. Some of the main 
results are as follows: (i) for a given RG equation the ex- 
plicit dependence on the regulator function disappears 
in the LPA, i.e. the results should not depend on the 
form of the regulator function if no further approxima- 
tions (e.g. truncation in powers of the field) are used, (ii) 
RG equations linearized around the trivial UV Gaussian 
fixed point (K = 0) provide the same UV scaling laws in 
various RG schemes, consequently, critical exponents of 
the Gaussian fixed point are scheme-independent in the 
LPA, (iii) non-trivial fixed points (like the Wilson-Fisher 
fixed point) and the critical exponents characterizing the 
scaling in their neighbourhood obtained in the LPA be- 
come scheme-dependent, i.e. the estimates of these phys- 
ical quantities depend on the particular choice of the RG 
scheme. Therefore, it is generally assumed that once ap- 
proximations are used non-trivial fixed points and their 
critical behavior become scheme-dependent. However, 
in this paper we show that one can find models which 
have important physical realizations and their physical 
parameters (fixed points, critical ratio) can be obtained 
scheme-independently even in the LPA. 

Our goal here is to consider the scheme-dependence 
of the low-energy behavior of various RG methods by 
considering the renormalization of the 2D periodic scalar 
field theory, i.e., the 2D SG model [SJ with possible inclu- 
sion of explicit symmetry breaking mass terms, i.e., the 
MSG [i| and the multi-component LSG model Q. The 
functional RG approaches are investigated in the LPA. 
Four types of RG methods are compared: the Wegner- 
Houghton (WH-RG) [1^, the Polchinski (P-RG) [H, 
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the Callan-Symanzik type (CS-R G) [3(]|| and the effective 
average action (EAA-RG) |3ll[3^.l33j| approaches. In the 
latter case we use two types of regulator functions: the 
optimized [6*1 and the power-like (quartic) [33j regulators. 
The phase structure of various SG type models obtained 
at quantum and classical level are also compared. 

As a rule, scheme-dependence is expected even in the 
exact (not truncated) RG flow at intermediate scales be- 
tween the UV and IR scales since various schemes re- 
alize the elimination of quantum fluctuations in a dif- 
ferent manner. However, the disappearance of scheme- 
dependence of the exact RG flow is expected in the deep 
IR limit if no approximations are involved. In the UV 
limit even the LPA of the RG flow is able to produce 
results which are independent of the particular choice 
of the renormalization scheme since various RG equa- 
tions linearized at the trivial Gaussian fixed point are 
the same in the LPA. Nevertheless, RG flows obtained 
in the LPA show spurious scheme-dependence at the 
non-trivial Wilson-Fisher fixed point. Therefore, it is a 
natural question to ask whether it is possible to obtain 
scheme-independent results at non-trivial fixed points. 
The answer can be affirmative if an additional constraint 
influences the RG flow. For example theories exhibiting 
spontaneous symmetry breaking have superuniversal ef- 
fective IR behavior due to the Maxwell-cut in the symme- 
try broken phase. Superuniversality represents so strong 
constraints on the RG flow that the scheme-dependence 
disappears if the effects of the truncation of the func- 
tional subspace are under control by using a sufficiently 
large functional subspace. In this paper, we will show 
that this general view holds for the MSG model where 
scheme-independent results are obtained for the critical 
ratio {u/AP)c being one of the typical IR characteris- 
tics of the Ising type phase transition of the MSG model. 
It will, however, also be shown that RG flow equations 
linearized at the UV Gaussian fixed point provide us 
scheme-dependent results for the critical ratio (u/M^)c 
of the MSG model which characterizes the phase struc- 
ture far beyond the validity of the UV scaling laws, in 
the deep IR region. 

Consequently, one expects that UV scaling laws cannot 
be used to determine the critical behavior of any non- 
trivial fixed point situated in the IR or in the crossover 
regions. At cross-over scales between the UV and IR 
scaling regions even "exact" RG equations (where LPA 
is the only approximation used) produce results being in- 
fiuenced by the choice of renormalization scheme. Never- 
theless, we shall show that the critical frequency = Stt 
at which the SG model and the layer-dependent criti- 
cal value (31{Nl) = %ttNl/{Nl - 1) at which the multi- 
component [Nl > 1) LSG models undergo a KTB type 
phase transition can be obtained exactly (independently 
of the invented RG scheme) even by the UV linearized 
RG fiow. This rather surprising result is due to the fol- 
lowing circumstances: (i) although UV linearized flow 
equations obtained in various RG schemes are different 
for the LSG model (where the linearization is performed 



in the periodic piece of the potential and not in the full 
potential), they provide us the same scheme-independent 
critical frequency, (ii) the KTB type phase transition 
of SG and LSG models is governed by the fundamental 
Fourier mode which has the same scaling law in the UV 
and IR regions, (iii) higher harmonics are generated by 
RG transformations and they have different UV and IR 
scalings but they do not influence the critical behavior. 

III. SINE-GORDON MODEL 

In this section we compare the applicability of the fre- 
quently used functional RG methods to two-dimensional 
one-component Euclidean scalar fields with periodic self- 
interaction. In Appendix we review the notation and 
the main properties of various RG schemes used by us 
which are the following: the Wegner-Houghton (WH- 
RG), the Polchinski (P-RG), the Callan-Symanzik type 
(CS-RG) and the effective average action (EAA-RG) ap- 
proaches. In the latter case we use the optimized and the 
quartic regulator functions. The WH-RG, the functional 
CS-RG, and the EAA-RG with power- like regulator with 
6=1 are equivalent methods in the LPA for d = 2 even 
if spinodal instability (SI) [13] occurs (see Appendix [X)) . 
Therefore, we shall compare the WH-RG, the EAA-RG 
with optimized and quartic regulators and the P-RG. It 
was shown [1 that the P-RG and the EAA-RG with 
optimized regulator can be transformed into each other 
by a suitable Legendre transformation (in LPA), how- 
ever, their singularity structures are different. We shall 
demonstrate this for the SG model since a SI, i.e. an IR 
singularity appears in the RG fiow in the weak coupling 
phase but the P-RG method is not able to indicate this. 
Therefore, we shall conclude that any of the discussed RG 
schemes are expected to be applicable to the investigation 
of the IR behavior of the MSG model except the P-RG 
method which is inappropriate for quantitative analysis 
in that case. 

For the generalized SG model characterized by the lo- 
cal potential 

N 

Vk{(t>) = Y. Un{k) cos {n/3(j,), (1) 
ji=i 

exhibiting periodicity in the internal space the dimen- 
sionless couplings are represented by the Fourier ampli- 
tudes Un{k) and the 'frequency' /9 is a scale-independent, 
dimensionless parameter in the LPA. If the higher har- 
monics are neglected, i.e. for TV = 1 the generalized SG 
model reduces to the well-known SG model 3]. We re- 
strict ourselves to the RG analysis of bare models with 
Ml (A) > (positive fugacity for 'charges' ±1 of the 
equivalent Coulomb gas). Another interesting generaliza- 
tion of the periodic model is the multi-component LSG 
model where the two-dimensional SG interaction terms 
are coupled by a particular mass matrix. This model re- 
ceives important applications in high-energy and in low- 
temperature physics and its KTB-type phase transition 
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has been discussed in [H, [Tl, [H, [H, [H, [111. In Sec- 
tion |V] we consider the scheme-dependence of the critical 
frequency of the multi-component LSG model. 

Let us first briefly summarize the results of the WH~ 
RG analysis of the generalized SG model obtained pre- 
viously in the LPA [13, [13, [H, [11] . The phase structure 
(in LPA) is sketched in Fig.[T]and has also been discussed 
in . The 0^ axis represents the line of Gaussian fixed 
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FIG. 1: The phase structure of the SG model obtained in LPA 
{0^ is scale-independent). The arrows indicate the direction of 
the RG flow. The non-trivial IR fixed points for < Stt are 
obtained by the WH-RG (full line) and by the P-RG (dashed 
line) methods. The line of IR fixed points is given by the 
analytic expression Mi(0) = 2//3^ for the WH-RG method for 
P'^ < Stt and ui(0) ~ 0.03(8% ~0y/'^ for the P-RG approach 
for ^8n. 

points, which are IR ones in the strong-coupling phase 
with > Stt and UV ones in the weak-coupling phase 
with < Stt. The two phases are separated by the verti- 
cal line at = Stt where the KTB type fll'l phase tran- 
sition takes place. In the figure only the single parameter 
axis ui is depicted. In order to get a more reliable picture 
of the phase diagram, one has to imagine a continuous 
sequence of non-intersecting infinite-dimensional hyper- 
surfaces which are intersected by the axis at {t„ = 0, 
as well as, by the line of the IR fixed points in the weak- 
coupling phase at some 'u„(0) ^ 0. Then the hypersur- 
face at = Stt separates the strong- and weak-coupling 
phases. In the UV limit {k < A) the linearized RG equa- 
tions provide the UV scaling laws 

u„(/c) = {i„(A) (^^^ " . (2) 

In the strong-coupling phase (/3^ > Stt) all the Fourier 
amplitudes are UV irrelevant whereas in the weak- 
coupling one (/3^ < Stt) the first few Fourier amplitudes 
become relevant, depending on the value of (3'^. In both 



the strong- and weak-coupling phases the dimensionful 
blocked potential tends to a constant effective potential 
for fc — > 0, being the only function which is simultane- 
ously convex and periodic, but there is a significant dif- 
ference in the IR scaling of the dimensionless couplings. 
The insertion of the ansatz ([T]) into the WH-RG equation 
Eq. (IA2p yields the RG flow equations 

{2 + kdk)nUn = — n^M„-f — ^s^„,s(2-l-fc9fc)Ms, (3) 

s— 1 

for the couplings u„ where An,s{k) = in — s)^u\n-s\ ~ 
{n-\- sYun+s- Eq- (|A2p . more precisely Eq. (|Aip are valid 
unless SI arises. In the strong-coupling phase > Stt 
no SI occurs, Eq. ([3]) holds at any scale and every Fourier 
amplitude is irrelevant. The IR scaling laws are given by 

/ A; \ 

Un(fc) = c„ f - j (4) 

with 77 = /3^/47r — 2 > and the constants c„ depending 
on 0^ and the single bare parameter ci = 'Ui(A) via the 
recursion relation (for n > 1) [3^ 

— TW\ ■ ^ ' 

n(2 + nr]- n^j:;^) 

Here the well-justified approximation An,s{k) « (n — 
s)^'S|„_s| has been used. The dimensionless blocked po- 
tential becomes flat, i.e. all couplings Un vanish in the 
IR hmit fc ^ 0. For ^2 <• 8^ ^j^g gj occurs in the RG 
flow when the propagator diverges, /c|j -I- Vl'^^{(j)) — 0. 
The scaling laws just above the scale kgi of the SI are 
given by (U, but now with rj < which modifies the 
scaling essentially. Namely, all the Fourier amplitudes 
Un become relevant at the scale ksi- Even more radical 
change of the IR scaling laws has been observed. Mak- 
ing use of the tree-level blocking relation (|A3|) one finds 
that the SI results in the building up of a condensate 
of increasing amplitude p for decreasing scale k < fcgi. 
In the weak coupling phase the dimensionless effective 
potential remains a non-vanishing periodic one (the line 
of non-trivial IR fixed points in Fig.[T]), graphically ob- 
tained by setting forth along the (p axis the section of the 
parabola 

2 °° C_L)"+i 1 
Vfe^o ^ E -^T- cosW) = - (6) 

n— 1 

with (f) G [— 7r//3, 7r//3] periodically. Each parabola sec- 
tion is the one of just the same parabola (|A4p . one would 
find as the non-trivial fixed point of the polynomial the- 
ory. The periodic dimensionless effective potential is the 
continuous, sectionally differentiable, periodic solution of 
Eq. (jA5[) . Let us note that Eq. ^ yields the fixed point 
equation 

2nUn = ^"^"i + E ^^■n^s'i.Us (7) 
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if no SI occurs in the RG flow. 

Using the same machinery in the framework of the P- 
RG, one obtains from Eq. (|A9|) the flow equations 

(2 + kdk)nun = —n^un ~ -0^ V sA„,^2ms- (8) 
47r z ^ — ' 

Inserting the ansatz ((4|) into Eq. ^ we get 

(2 + nry)nc„fc"'' = ^n^^™^"'' - ^ V sA„,,2c,r'' 

s— 1 

which gives 77 = 0^ /Att — 2 again. For > Stt ansatz 
Q works with the neglection of the term proportional 
to {t„_|_s(fc) in the expression of A„^s(fc) similarly to the 
WH-RG and yields the recursion relation 

with ci = ui(A), c„ = {t"(A)i?„ for 71 > 1, i?i = 1 and 
the constants i?„ for n > 1, satisfying the recursion re- 
lation ([9]), become independent of the bare couplings. In 
Fig. [2] the flow of the first few couplings can be followed. 
All dimensionless couplings tend to zero in the IR limit 
A; — » again. 




k/A 

FIG. 2: The P-RG flow of the first few dimensionless cou- 
pfings of the generalized SG model for = 127r. 

For 0^ < Stt the numerical solution of the flow equa- 
tions ([5]) for the dimensionless couplings Un{k) follows 
the power law behavior of Eq. (|4]), now with 77 < 0, in 
a wide range above a certain scale fee, but the couplings 
Un{k) go to constant values in the deep IR region, as it 
is demonstrated in Fig. [31 Let us note that the inverse 
propagator vanishes around the scale ~ kc which signals 
the appearance of the SI (i.e. kc ^ fcsi)- However, the P- 
RG flow equations do not exhibit any singularity at that 
scale kc- Eq. ([5]) provides a scale independent, IR fixed 




k/A 

FIG. 3: The flow of the first few dimensionless couplings of 
the SG model for 13^ — Att obtained by the P-RG method. 

point solution now satisfying Eq. ([7]) . Our numerical re- 
sults show in Fig.[3]for close to but below Stt < Stt) 
that the /3-dependence of Un{0) can be factored out as 

<(0) = (8^-/32)»/2a„, (10) 

where the numbers a„ are independent of any parameter 
of the model. Therefore, the non-vanishing, universal 
dimensionless IR effective potential obtained by the P- 
RG method reminds one on the similar result of the WH- 
RG analysis, but the parabolic shape of the dimensionless 
effective potential in its periods cannot be recovered. As 
a consequence, the line of IR fixed points of the SG model 
in the weak coupling phase has been modified, see the 
dashed line in Fig.[TJ 

As to the ambiguity of the decision based on the nu- 
merics whether the SI does or does not occur during the 
fiow, it is in order to make here an important remark. 
At the first glance one expects that the accuracy of such 
a decision will be increased by increasing the number of 
Fourier-modes, that of the couplings Un taken into ac- 
count. With decreasing scale k the dimensionful peri- 
odic potential becomes rather flat in each period. Ac- 
cording to our numerical experience it approaches rather 
smoothly the potential (O for ~ An and then the oc- 
currence of the SI can be detected without any doubt. 
However, the matter of things becomes much worse for 
< Stt — Z?'^ ^ Stt when there occur numerical instabili- 
ties if one reconstructs the potential. In such a case one 
cannot decide unambiguously with the numerical method 
based on the Fourier expansion of the potential whether 
the SI docs occur indeed. Solving the RG equation de- 
rived in the LPA for the potential without using any fur- 
ther expansion seems to be a reliable way to settle this 
point, but it lies out of the scope of the present paper. 

As to the next, let us turn now to the discussion of 
the EAA~RG flow for the generalized SG model (P). Us- 
ing the optimized regulator, and deriving Eq. (jAlSp with 
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FIG. 4: The IR values Un{0) of the dimensionless couphngs 
show power law behavior in the vicinity of the critical value 
of The markers are obtained by solving the non-linear 
system of equations in Eq. (O, and the lines are plotted by a fit 
giving 0.5, 1.0, 1.5 for the slopes, respectively for n — 1, 2, 3. 



respect to the field and multiplying its both sides by 
(1 + V^")^; one obtains tlie evolution equations 

(2 + kdk)nn^ = £^n^< + E^i P^n,p(2 + kdk)up 
-3/3' E^i Ep=lPAn,,A,j2 + kdk)u, (11) 

for the Fourier amplitudes Un- For > Stt the IR scaling 
laws are given again by Eq. ^ with -q — (3"^ / (47r) — 2 > 
and the recursion relation for c„'s 

_ Sp=i (2 + pr])p{n - pfcn-pCp 

n(2 + nr] - ri^j:;^) 

For 0^ < Stt the numerical solution of the system ([TT|) of 
the coupled flow equations exhibits the following features. 
If one takes into account a sufficiently large number of 
Fourier modes, one finds that there exists a non- vanishing 
scale fcgi at which the inverse propagator vanishes that 
signals the presence of the SI. The numerical problem of 
deciding whether the SI occurs or does not, remains just 
the same as in the WH-RG framework. For in the 
vicinity of the critical value Stt and restricting oneself to 
the first few Fourier amplitudes, one finds that the inverse 
propagator does not vanish in the IR region. In this 
case the IR effective potential is similar to that obtained 
by the P-RG, with the scaHng property pU]) (but with 
different numbers a„). 

Let us now discuss the flow in the framework of the 
EAA-RG with quartic regulator. For f3'^ > Stt the SI 



does not occur. On the basis of the result of the WH-RG 
analysis, for < Stt one would expect the occurring of 
the SI, i.e. vanishing of 1 + ^V^' {(f)) at some scale fcsi 7^ 0. 
For the single mode potential Vk{(j>) = u(k) cos(/30) this 
happens for (/> = if 1 — ^u(fcsi)/3'^ = and then 1 — 
iy"^"^ (tt//?) — holds as well. Therefore, the momentum 
scale kj^ where the quartic regulator RG is non-analytic 
and the scale of the SI fc_ = k^i coincide for the periodic 
model. The decision based on numerics on the existence 
of the SI suffers from the same problems as for the other 
RG schemes, and the IR scaling laws are also expected 
to be similar to those obtained on the basis of Eqs. ([3]), 
(III), and (fTTj) . Since the SI occurs when the propagator 
in the right hand side of Eq. (|A13p develops a pole, and 
periodicity should not be violated, the IR fixed point 
potential should be obtained by setting forth periodically 
the parabola section in [— 7r//3, 7r//3] given by Eq. (|A14p 
and Eq. (|XT5)) . 

In conclusion, the IR scaling laws determined by vari- 
ous RG methods are qualitatively the same in the strong 
coupling phase of the generalized SG model. In the 
framework of WH-RG and the EAA-RG with the power- 
law regulator, the SI has been treated explicitly in the 
weak coupling phase which provides a reliable determina- 
tion of the IR scaling laws even in that phase. In this case 
the functional form of the low-energy effective potential 
is found to be the same but its exact value is scheme- 
dependent (see e.g. Eq. (|A14[) and Eq. (|A15p ). However, 
the potential can always be rescaled by a constant which 
leaves the physical results unchanged, consequently, only 
the functional form of the IR effective potential is of phys- 
ical significance. 

IV. MASSIVE SINE-GORDON MODEL 

The MSG model for the number of dimensions d = 2, 
characterized by the dimensionless bare potential 

Vk=K{<t>) = \m14>'' + «(A) cos(/30) (13) 

exhibits two phases. For 0^ > Stt only the phase with ex- 
plicitly broken periodicity is present, and that phase ex- 
tends to < Stt if the bare dimensionless ratio i2(A)/M^ 
is smaller than a critical upper bound ['u(A)/M^]c de- 
pending on the parameters and M^. For < Stt and 
u[h)/M\ > [u(A)/M|]c the periodicity is spontaneously 
broken, and the RG trajectories in that phase merge into 
a single trajectory in the deep IR region which is charac- 
terized by the unique ratio u{k 0)/Af|_^Q = [u/M'^]c = 
[u/M'^]c depending on 0^ only. Our main goal is to dis- 
cuss whether the determination of that unique ratio does 
depend on the choice of the previously discussed renor- 
malization schemes. Let us note that the Z2 symmetry 
(j) ^ —(j> can be used to distinguish the phases of the 
MSG model, it suffers a spontaneous breakdown in the 
phase with spontaneously broken periodicity when a con- 
densate appears. 



7 



The phase structure of the MSG model has been dis- 
cussed in the framework of the WH-RG in the LPA in 
[l3| making use of the more general ansatz 



1 ^ 



(14) 



3. Inserting the ansatz (|T4|) into the WH-RG equation 
Eq. (|A2p yields the RG flow equations 



M2 



13' 



(1 + —){2 + kdk)nun = ^n'^Un + (18) 



for the local potential. (The parameter u in the ex- 
pression corresponds to ui(A) and we are looking 
for the unique value of ui{k — > 0)/Af^^Q = [u/M'^]c 
in the deep IR region for the phase with spontaneously 
broken periodicity.) The periodic piece of the poten- 
tial (fT4|) possesses the discrete symmetry under the shift 
(j) ^ (j) + 2tt/P of the field variable. This is the sym- 
metry the generalized SG model with the potential ([1]) 
exhibits, and - as discussed in the previous section - in 
the weak couphng phase (/3^ < Stt) of the generahzed SG 
model a condensate appears in the IR limit, due to which 
the periodicity is broken spontaneously. Due to the ex- 
plicit mass term the MSG model offers the opportunity 
to investigate the interplay between the spontaneous and 
explicit breaking of periodicity. The following features of 
the phase structure have been determined in the frame- 
work of the WH-RG method in the LPA in [13, HI: 

1. The dimensionful mass remains constant dur- 
ing the blocking which provides the trivial scaling 

= k~'^A'P for the dimensionless mass. In the 
LPA (when no wave-function renormalization is in- 
corporated) also the parameter 0^ is constant dur- 
ing the flow. 

2. The linearized WH-RG flow equation 



{2 + kdk)ui{k)^^ui{k)- ''^ 



47r 'k^ + A'P 



(15) 



for the couphngs Un{k). For fc^ ^ the RG 
flow is close to that of the SG model. For strong 
coupling > Stt no SI occurs during the flow, 
the inequality 1 + V^' > holds for all scales k. 
The numerical solution of Eq. ^TE\i provides scaling 
laws for M < fc ^ A not differing significantly 
from those of the SG model, and below the scale 
k ~ 0{M) the trivial scaling Un{k) ~ k~'^ occurs. 
The parameter region with 0^ > Stt belongs to the 
phase with explicit breaking of periodicity. 

4. For weak coupling 0^ < Stt two phases appear. 

On the one hand, the RG flow on the trajectories 
started at /J^ < Stt and mi(A)/M^ > [ui{A)/Ml]c 
develops SI at some finite scale fcsi > M, so that 
one has to evaluate the flow using the tree-level 
blocking relation IjASp . In this case the parameters 
Un{k) become superuniversal (i.e. independent of 
the bare parameters) for k < kgi signaling the pres- 
ence of the condensate and the spontaneous break- 
ing of periodicity, similarly to the behaviour of the 
weak coupling phase of the SG model. It should 
be noticed that the detection of occurring SI in the 
evaluation of the trajectories suffers the same am- 
biguity for close to but smaller than Stt as in the 
case of the generalized SG model. For k < ksi the 
tree-level evolution yields the periodic piece of the 
potential with parabola sections. 



for the coupling ui{k) of the fundamental mode of 
the local potential provides the UV scaling law 



Ul 



(fc)=ui(A) - 



A / V A2 + M2 



(16) 



The coupling ui{k) is relevant in the UV scaling 
region, so that the scaling law (|16p looses its valid- 
ity for IR scales, irrespectively of 0. Nevertheless, 
it was observed numerically that Eq. (fT6|) provides 
a rather good description of the RG flow for scales 

k^ M for > Stt. The critical ratio 



Ml (A) 

Ml 



1 



Ml 



2Ml 



(17) 



for the bare parameters of the MSG model has also 
been discussed in the framework of the linearized 
RG (see Eq.(9) of 0). 



^ n=l 

= -]^{k' + M')(t? 



cos{nf3(j)) 



(19) 



for <j) e [— 7r//3, 7r//3]. Consequently, the first di- 
mensionful Fourier amplitude tends to the constant 
ui{k ^ 0) = ^M^, i.e. the IR value of the critical 
ratio is [u/A/^J^ — 2/0^. The parabolic potential 
(|19|) represents the non-trivial solution of the differ- 



ential equation F + Vfc"(0) = k^ + M^ + U'^((j)) = 0. 
The effective potential Vk^o — Uk^o + ^M'^cj? = 
— ^k'^cj)'^ is superuniversal, i.e. independent of the 
bare parameters u„(A). 

On the other hand, the phase with explicit breaking 
of periodicity extends to the region with (3"^ < Stt 
bounded from above by [ui{A)/Ml]c. In this phase 
no SI occurs along the RG trajectories and below 
the scale k = M the Fourier amplitudes start to 
scale as M„(fc) k~^. In this case the effective 
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potential Uk^o turns out to depend on the single 
bare parameter ui(A). 

The phase diagram of the MSG model for — 4tt is 
depicted in Fig. [5] where one can see that the RG trajec- 
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FIG. 5: Phase structure of the MSG model in LPA for (fi = 4tt 
obtained by the full WH-RG method. The arrows indicate 
the direction of the flow. The numerical solution of the ex- 
act RG equation shows that the separatrix obtained by the 
linearized flow (I21|l (wide dashed line) is modified by higher 
order corrections. 

tories in the phase with spontaneously broken periodicity 
merge into a single trajectory in the deep IR limit (full 
lines) . It is worthwhile noticing that the KTB phase tran- 
sition exhibited by the massless SG model disappears in 
the MSG model due to the presence of the explicit mass 
term, as demonstrated in [13, 123 using WH-RG and CS- 
RG methods. 

Let us consider the critical RG trajectory on the ui — 
plane which separates the phases of the MSG model. 
At linearized level this is constructed from Eq. as 



1 + 



A, 



1 



Ml 



(20) 



where = M| = k "^A'P is the running coupling and 
Ml = A~^M^ is constant. By inserting the critical ratio 
P?)) obtained for the initial bare values into Eq. (|^ one 
finds 



UlciM') = —M 



1 + M^ 
2M2 



At the mass-scale it gives uic{M'^ = 1) = 2//?^ 
(see the wide dashed line in Fig. [5]). From Eq. 
can read off the IR value of the critical ratio 



M2 



= lim 



Af2 



2 



(21) 

w 0.159 
2T]) one 



(22) 



which is smaller than the superuniversal ratio [u/M'^]c = 
2//32 obtained by the exact WH-RG method (for /J^ ^ Q) 
and coincides with it only in the limit 0^ — > 0, see the 
dashed line in Fig. [SI Indeed, the numerical solution of 
the exact WH-RG equation (|18p shows that one finds 
RG trajectories (see. Fig. [5]) which lie above the critical 
one obtained in the framework of the linearized WH-RG. 

In Fig.[B]we plotted the superuniversal ratio [u/M'^]^ = 
ui{k 0)/M2 against the parameter f]^ obtained in 
the framework of the WH-RG. According to our numer- 
ical results the effective potential implies [m/M^]^ = 
2//52. According to lattice results [19| the phase transi- 
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FIG. 6: The critical ratio which separates the phases of the 



MSG model obtained at classical level {u/M 



1/(3-' 



compared to that of determined by WH-RG and EAA-RG 
methods. The full circle represents the results (u/M^)c = 
0.156 — 0.168 calculated by lattice methods [l9l . [20l . [2ll ] at 



4n. The solid line stands for {u/M^)^ = 2//3'^ which 



is obtained by the full WH-RG (-1- tree-level RG) method. 
(The full EAA-RG (-1- tree-level RG) method provides the 
same critical ratio (solid line). The dashed line is determined 
by the linearized WH-RG equation and it is given by (|22|l . 
The dashed-dotted fine is obtained by the linearized optimized 
EAA-RG approach, see Ea. (|27|l which gives better result then 
the linearized WH-RG. For 0^ = 4-k it is shown that only 
the usage of the full RG equation gives reliable result for the 
critical ratio i.e. {u/M'^)c = 0.159 which coincides to the 
results of the lattice calculations. 

tion of the MSG model is assumed to belong to the same 
universality class as the two dimensional Ising model 
(with the critical exponents i^ = l,/3=l/8 [1^). Since 
the MSG model for (3'^ = Air represents the bosonized 
version of QED2, the ratio [u/M^lc = l/(27r) w 0.1591 
obtained by the WH-RG method can also be used to 
determine the phase transition point of QED2 given as 
[m/g]c = 2^e~'^[u/M^]c ~ 0.3168 where m is the 
fermion mass, g is the coupling between the fermionic 
and the gauge field and 7 w 0.5774 is the Euler's con- 
stant. It is shown in Fig. [5] that our WH-RG result is in 
good agreement with the result [m/g]c = 0.31-0.33 ob- 
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tained by density matrix and lattice methods [l9|, |20|, |2l| . 

It is illustrative to compare the critical ratio [m/M^]c = 
2//32 obtained by the WH-RG method with the esti- 
mation based on the classical potential p3p by solv- 
ing V;i{(l) = 0) = which gives [u/M\ = l//?^ i.e., 
[m/g\c — 2^e^'^ 1 0^ k, 0.1584, a value being one-half 
of the one obtained by taking the quantum fluctuations 
into account. This means that the slope of the critical 
RG trajectory which separates the phases of the MSG 
model is steeper when the effect of the quantum fluc- 
tuations is taken into account, i.e. the latter enlarges 
the parameter region with explicitly broken symmetry, 
see Fig. [71 In general, the effect of quantum fluctuations 



classical lineaiized-RG 



full-RG 

V 



explicitly broken periodicity 



spontaneously 
broken periodicity 



linearized evolution equation 

/32 



(2 + kdk)ui 



M2)2 



(23) 



for the fundamental Fourier amplitude when the opti- 
mized regulator is made of use. Here and /J^ are 
again scale-independent parameters and the analytic so- 
lution. 



ui{k) = ui{A) ( - 



exp 



M 



,t2 



A2 + Af 2 

M2 



fc2 -f M2 A2 + A/2 



(24) 



provides qualitatively the same UV scaling law as that of 
Eq. pB)) and indicates the absence of the KTB-type phase 
transition. Let us determine the critical RG trajectory 
on the ui ~ AP plane at linearized level. 



FIG. 7; Regions of explicitly and spontaneously broken pe- 
riodicity of the MSG model obtained by various approxima- 
tions. The critical value of the fundamental Fourier amplitude 
ui{ki-R.) is calculated by classical, linearized- RG and fuU-RG 
methods at some IR momentum scale fciR, << M << A. The 
effect of quantum fluctuations leads to the shrinking of the 
region of spontaneous breakdown of symmetry. 



leads to the shrinking of the region of spontaneous break- 
down of symmetry. For example, the polynomial model 
defined by Eq. (jBip exhibits a double-well bare potential 
for the initial conditions with 52(A) < chosen in Ap- 
pendix |b1 but the running mass 52 (fc) changes its sign 
due to quantum effects (see Fig. [8]) and the minimum of 
the effective potential remains at = 0. Let us note that 
the classical analysis is not able to distinguish between 
the weak {(3'^ < Stt) and the strong {(3"^ > Stt) couphng 
regimes, it determines the critical ratio [u/M^]^ = 1/(3"^ 
independently of the actual choice of the frequency (3'^. 
Only the detailed RG study of the MSG model can show 
that for (3'^ > Stt the SI, i.e. the condensate does not 
appear, hence the critical value cannot be extended for 
/32 > Stt, see the vertical line in Fig.[6l 

In Appendix |B1 it was demonstrated that beyond the 
mass scale in the deep IR limit the RG flow of the poly- 
nomial scalar field theory (in the phase with unbroken Z2 
symmetry) determined by the P-RG differs of the usual 
one obtained by other RG methods. It was also shown in 
Section IIIII that the P-RG approach is not able to signal 
the appearance of SI for the weak coupling phase of the 
SG model. Consequently, the P-RG is inappropriate to 
investigate the IR scaling of the MSG model quantita- 
tively, hence, we do not discuss the RG flow of the MSG 
model in the framework of the P-RG method. 



Let us turn to the discussion of the MSG model (|M1) in 
the framework of the EAA-RG. Considerations similar 
to those made in the case of the WH-RG lead to the 



ui,(M2) = ui,(A) 



M2 



I + Af2 



I 



exp 



Af2 



f Af2^ 
Ml 



1 + A/2 I 



Ml, 



(25) 



The critical initial bare value uic{A) is fixed by the condi- 
tion uic(A/2 = I) = 2//32 similarly to the WH-RG case, 
then one finds 



exp 




(26) 



From 
ratio 



one can read off the IR value of the critical 



A/2 



= lim 



Ulc 
^°°A/2 



2 



eSTT 



(27) 



which is smaller then the superuniversal ratio [u/M'^]c = 
27/32 obtained by the exact WH-RG method but gives 
a better result then Eq. ([22]) obtained by the linearized 
WH-RG equation, see the dashed-dotted line in Fig.[Sl 
Indeed, the optimized regulator has been constructed to 
achieve the best convergence of the truncated RG equa- 
tion, consequently, between various RG equations consid- 
ered at the same order of the truncation, the optimized 
one gives the closest result to the exact one [6] . 

However, neither the linearized WH-RG nor the lin- 
earized optimized EAA-RG enables one to map the phase 
structure of the MSG model in a reliable manner, i.e. to 
determine the exact critical ratio [u/A/2]^ which can only 
be obtained by solving the full RG equation in the IR 
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limit k ^ 0. As shown in the framework of the WH-RG 
method, if SI arises, its appropriate treatment is neces- 
sary to determine the effective potential beyond the scale 
of SI and the critical ratio [u/M'^]c- The flow equation 
(|A18|) has a pole at Uj! + IvP = — fc^, consequently, the 
IR effective potential beyond the scale fcsi is found to be 
identical with that of Eq. (|X4)) obtained by the WH-RG 
method. Therefore, the critical ratio [u/M^]c coincides 
with that determined by the WH-RG approach (see the 
solid line in Fig.E]). 

Finally, in the framework of the EAA-RG with the 
quartic regulator the linearization of Eq. (jA17[) results in 
the flow equation 



{2 + kdk)ui = ^^iT^ 



4 



(28) 



l + Mfe(arctg/i2-| 
1 A ^„ ( 



for \1VP < 
for > P 



with nl = iMV^|/c4 _ iM4| and the scale- 
independent parameters and Eq. P5|) provides 
again an UV scaling law which is very much like those 
of Eq. m]) and Eq. ([M]). In order to find the vahdity 
range of the mass-corrected UV scaling law and to map 
the phase structure, one has to solve Eq. (|A17[) numeri- 
cally. For sufficiently large initial values of mi(A) an IR 
singularity appears in the RG flow at some scale fcsi, sim- 
ilarly to what happened in the framework of the WH-RG 
and the EAA-RG with the optimized regulator. One has 
to use Eq. (|XT6)) . i.e. + [/^' = ~C{b)P to treat the 
SI explicitly and arrives at Uk = -\[C{b)k'^ + M^](j)'^ 
and the IR effective potential (|A15|) with C(2) — 2, i.e. 
Vk^o = — fc^0^. Since Uk=o = — ^Af^^^ is independent 
of the scheme-dependent constant C{b), the critical ratio 
turns out to be independent of the scheme as well. Fur- 
ther on the potential Uk=o and, consequently the critical 
ratio [mi/M^]c = 2//3^ are just the same as those found 
by the WH-RG method. 

One can conclude, on the one hand that the WH-RG, 
the CS-RG (being equivalent now with the WH-RG), 
and the EAA-RG methods enable one to determine the 
same phase structure and critical ratio [m/M^]c for the 
MSG model. On the other hand the P-RG method is 
found to be inappropriate to follow the RG trajecto- 
ries beyond the mass scale and to perform a quantitative 
analysis of the IR behavior of the MSG model. It was 
also shown that the truncated RG equations produce a 
scheme-dependent critical ratio but any exact RG (LPA 
is the only approximation used) gives the same (scheme- 
independent) result which coincides with the results of 
density matrix and lattice methods. 



V. LAYERED SINE-GORDON MODEL 

In this section we show that the critical frequency 
which separates the phases of an SG-type model which 



undergoes a Kosterlitz-Thouless-Berezinskii type phase 
transition similar to the one for the two-dimensional 
SG model can be obtained exactly (and scheme- 
independently) by the linearized RG flow. As an example 
we consider the renormalization scheme-dependence of 
the LSG model, i.e. the multi-component scalar field the- 
ory where the two-dimensional periodic interaction terms 
are coupled by an explicit mass matrix 



Vkif) = If Kit + "(^) E cos(/30„) 



(29) 



n=l 



where is the number of the coupled fields (i.e. number 
of "layers") with the 0{Nl) multiplet (j) = (0i, ■ ■ • , 4>Nl)- 
The mass-matrix describes the interaction between the 
fields and is chosen here to be of the form 



(30) 



where Gk is the strength of the inter-field interactions. 

The WH-RG equation in LPA for the multi- 
component LSG model presented in our previous pub- 
lications fH, m [13, in reads as 




(31) 



where Vj^-' (</>) denotes the second derivatives of the po- 
tential with respect to Inserting the ansatz (|29|) 
into the RG equation (pij) the right-hand side becomes 
periodic, while the left-hand side contains both periodic 
and non-periodic parts. The non-periodic part contains 
only mass terms, so that we obtain a trivial tree-level RG 
flow equation for the dimensionless mass matrix 



{2 + kdk)K^ = o, 



(32) 



which provides the trivial scaling Gk — k^^G, where 
the dimensionful inter-field coupling G remains constant 
during the blocking. We recall that in LPA there is no 
wave-function renormalization, thus the parameter f3 also 
remains constant during the blocking. Although the so- 
lution of Eq. (|3ip can only be obtained numerically, how- 
ever, analytical results are also available using an approx- 
imation of Eq. pT|) similarly to the MSG model. This 
is achieved by linearizing the WH-RG equation in the 
periodic piece of the blocked potential (not in the full 
potential), 

{2 + kdk)Uk{c^i,...,cj)Nr.)~-^ (33) 



An C 



where Uk{4'i, 4>Nl) = u{k) Yl,n=i cos(/3 0„) and C and 
Fi{Uk) stand for the constant and linear pieces of the 
determinant 

det[J„- + V^'] « C + Fi{Uk) + 0{Ul). (34) 
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The mass-corrected linearized WH-RG for the coupled 
periodic model (|29p reads as 

where G and 0^ are scale-independent parameters and 
the solution can be obtained analytically 



u{k) = u{A) ( - 



where u{A) is the initial value for the Fourier amplitude 
at the UV cutoff A. The critical frequency which sepa- 
rates the two phases of the model can be read out directly 
as 



SttNl 
Nl~1' 



(37) 



For Nl = 1 the coupled model ([29]) reduces to the mas- 
sive 2D SG model with 0^ = oo which indicates the ab- 
sence of the KTB phase transition. For Nl — > oo the cou- 
pled SG model behaves like a massless 2D-SG model with 
the critical frequency — Stt. Our goal is to consider 
whether the determination of the critical frequency ((37|) 
does depend on the choice of renormalization scheme. 

Since the P-RG method fails to determine the correct 
IR behavior of the one-component MSG model, here we 
do not apply the P-RG method to study the properties 
of the multi-component LSG model. Therefore, let us 
turn directly to the RG analysis of the LSG model (|29|) 
by the EAA-RG method with the optimized regulator. 
Then the RG flow equation for the local potential of the 
A^i-component scalar field reads as 



(2 + kdk)Vk - -^Tr 
47r 



(38) 



where V^'' — d(p-d^-Vk. Using the same machinery as in 
case of the WH-RG, the mass-corrected linearized form 
of the optimized regulator RG for the LSG model (129]) 
reads as 

(3^ ^k^ + 2{Nl - l)Gfc2 + Nl{Nl - l)G^ 



(2 -f kdk)u = — u 
Ait 



(fc2 NlGY 



(39) 

where G and 0^ are scale-independent parameters and 
the solution can be obtained analytically. 



i(fc) = w(A) 



exp 



k 
A 

Stt 



NlG 



G 



A2 + NlG 
G 



fc2 + NlG A2 + NlG 



(40) 



where m(A) is the initial value for the Fourier amplitude 
at the UV cutoff A. This gives the same A''i,-dependent 
critical frequency ([57)1 as that obtained by the WH-RG 
method. 



In general, any RG method which is applicable to pre- 
dict the IR behavior for the one-component massless and 
massive SG models, hke the WH-RG, the EAA-RG with 
either the optimized or the quartic regulators, and the 
functional CS-RG methods, is assumed to produce the 
same A^L-dependent critical frequency ([37]) for the LSG 
model (l^ni) . This can be understood by using a suitable 
0{Nl) rotated form of the coupled LSG model where 
the mass matrix is diagonal. Since the rotation leaves the 
phase structure unchanged, therefore, instead of the orig- 
inal model one can investigate the rotated one where the 
single massive 2D-SG field can be considered perturba- 
tively which results in an effective SG-type model. This 
strategy can be applied for the multi-component cou- 
pled LSG model with arbitrary number of components 
{Nl > 1), hence, in the lowest order of the perturbative 
treatment the corresponding effective theory is always an 
SG-type model and, consequently, the various linearized 
RG equations produce the same critical frequency. 



VI. SUMMARY 

The renormalization of sine-Gordon (SG) type pe- 
riodic scalar field models with explicit mass terms, 
i.e., the massive sine-Gordon (MSG) and the multi- 
component layered sine-Gordon (LSG) models have been 
investigated by various functional renormalization group 
(RG) methods using the local potential approximation 
(LPA). Our aim was to compare the Wegner-Houghton, 
the Polchinski, the functional Callan-Symanzik methods 
and the effective average action RG method with var- 
ious regulator functions and to investigate the scheme- 
dependence of the low-energy behavior of SG type mod- 
els. In particular, our goal was to consider under 
which conditions is it possible to determine scheme- 
independently three physical parameters, the critical fre- 
quency of the SG theory, the critical ratio u/M"^ of 
the MSG model and the layer-dependent critical value 
(31{Nl) of the LSG model. 

Even if the effects of the truncation of the functional 
subspace in which the effective potential is sought for 
are under control, the RG flow depends on the particular 
choice of the renormalization scheme due to the different 
ways the quantum fluctuations are eliminated. Never- 
theless no scheme-dependence is expected if the UV scal- 
ing laws or the (exact) IR physics are considered. UV 
scaling laws are invoked from the linearized flow equa- 
tions being insensitive to the choice of the renormaliza- 
tion scheme. The IR physics is obtained by integrating 
out all the quantum fluctuations and the result does not 
depend on the manner how this happened if no approx- 
imations (e.g. a truncation of the functional subspace) 
were used. 

We demonstrated that the RG flows of the Wegner- 
Houghton, the functional Callan-Symanzik and the effec- 
tive average action RG methods are similar for polyno- 
mial models (in the phase with unbroken Z2 symmetry). 
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for the generalized SG model (in its both phases) as well 
as for the MSG model. Those methods enable one to de- 
termine the IR physics, the phase structure of the MSG 
model in a reliable manner, in good agreement with den- 
sity matrix and lattice results. It was also shown that 
the Polchinski RG method is inappropriate to determine 
the RG flow beyond the mass scale in the deep IR limit. 
While this is not a drawback in case of polynomial mod- 
els and the generalized SG model (with vanishing mass) 
it disables one to determine the phase structure of the 
MSG model quantitatively. 

The phase structure of the MSG model (which is the 
Bose form of QED2) obtained by classical and quantum 
analysis were also compared. It is known that QED2 has 
two phases, the strong g » m and the weak g << m 
coupling phase, where m is the fermion mass and g is 
the coupling between the gauge and the fermion field. 
The critical ratio [m/g]c which separates the phases of 
QED2 obtained at the classical level was found to be 
two times smaller than that determined by RG meth- 
ods. Above (below) the critical ratio, the periodicity of 
the bosonized model (i.e. for the MSG model) is bro- 
ken spontaneously (explicitly). Therefore, the region of 
spontaneous symmetry breaking in the parameter space 
is reduced at quantum level as compared to the classical 
one which is in agreement with the general assumption, 
that the effect of quantum fluctuations always shrinks 
the region of spontaneous symmetry breaking. 

On the one hand, we showed that the critical frequency 
which separates the phases of an SG-type model which 
undergoes a Kosterlitz-Thouless-Berezinskii (KTB) type 
phase transition similar to the two-dimensional SG model 
can be obtained exactly (and scheme independently) by 
the RG flow linearized around the UV Gaussian fixed 
point. This is the consequence of the extension of the 
UV scaling region down to the vicinity of the crossover 
region separating the UV and IR scaling regions. As 
examples we considered the renormalization of the one- 
component two-dimensional SG theory and the multi- 
component LSG model where the two-dimensional SG 
fields of the various layers are coupled by an explicit mass 
matrix On the other hand, it was also shown that 

the linearized RG equations produce a scheme-dependent 
critical ratio [u/M'^]c for the MSG model (i.e., [m/g]c 
for QED2) and only the "exact" RG (where LPA was 
the only approximation used) gives the same (scheme- 
independent) result which coincides with the critical ra- 
tio determined by density matrix and lattice methods. 
Therefore, this demonstrates that a KTB-type phase 
transition of an SG-type model is a crossover between 
the UV and IR scaling regimes where the linearization 
around the Gaussian UV fixed point gives reliable results 
but the Ising-type phase transition of the MSG model is 
found to be an IR one where one has to solve the full 
RG equation in order to determine the exact transition 
point. 

Finally, let us conclude with two general results which 
are not restricted to the particular models investigated. 



We have shown that in some cases it is possible to ex- 
tend the validity of the UV linearized RG flow down to 
the vicinity of the crossover region separating the UV and 
IR scales which can produces physical parameters (such 
as the critical frequency of the SG and the layer num- 
ber dependent critical value 01{Nl) of the LSG models) 
independently of the particular choice of the renormaliza- 
tion scheme. This is opposed to the general assumption 
namely that once approximations are used non-trivial 
fixed points and their critical behavior should be scheme- 
dependent. It was also shown that in case of spontaneous 
symmetry breaking the Maxwell construction represents 
a strong constraint on the RG flow which results in a su- 
peruniversal IR behavior and as a consequence, scheme- 
independent IR results (such as the critical ratio [u/M'^]c 
of the MSG model) can be obtained even if the local po- 
tential approximation is used. This receives important 
application in any case where the low energy effective the- 
ory can only be determined by non-perturbative methods 
like functional RG approaches and a spontaneous sym- 
metry breaking infuences the low energy behavior. Our 
results indicate that scheme-independent results can be 
obtained even in the LPA. 
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APPENDIX A: RENORMALIZATION GROUP 
METHODS FOR d = 2 IN THE LOCAL 
POTENTIAL APPROXIMATION 

In this appendix we settle our notations and remind 
the reader on some well-known features of the frequently 
used renormalization schemes in the local potential ap- 
proximation (LPA) for two-dimensional Euclidean one- 
component scalar field theories. In principle the vari- 
ous renormalization schemes are constructed in such a 
manner that the RG flow starts at the bare action and 
provides the effective action in the IR limit, so that the 
physical predictions (e.g. the critical exponents) are inde- 
pendent of the renormalization scheme particularly used. 
Nevertheless, scheme-dependence may appear in the RG 
flow at intermediate scales due to the different manners 
the quantum fluctuations are eliminated in the various 
RG approaches. But even the physical predictions at 
the IR scales may depend on the used renormalization 
scheme if additional approximations are involved like the 
improperly strong reduction of the functional subspace 
in which the local potential is sought for or the lineariza- 
tion of the RG flow equations at the Gaussian fixed point. 
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Therefore, it is of relevance to clarify how far the results 
obtained are independent of the particular choice of the 
renormalization scheme used. 

The use of the sharp cutoff, i.e. the Wegner-Houghton 
(WH) RG approach [2811 makes the blocking transforma- 
tion transparent and simple since the modes to be elim- 
inated are well-defined. The price of this clarity is the 
incompatibility with the gradient expansion. One possi- 
ble solution for this problem could be the usage of the 
smooth momentum cutoff, where the higher frequency 
modes of the field are suppressed partially, but not elim- 
inated. This can be realized by Polchinski's [2§| method 
(P~RG) for the bare action and by the effective average 
action (EAA) RG approach [sil, [H, HI] with various 
types of regulator functions. As a rule, the solution of 
the truncated RG flow depends on the particular choice of 
the regulator. Another treatment of the RG which han- 
dles the effective action obtained by a suitable Legendre- 
transformation is represented by the functional Callan- 
Symanzik (CS-RG) or also called the internal space RG 
approach [30j , where the quantum fluctuations are sepa- 
rated according to their amplitudes instead of their fre- 
quencies or length scales. In this case the blocking pro- 
cedure is performed in the space of the field variable and 
not in the space-time, which is understood as the external 
space. 

One of the advantages of the Wilsonian RG method 
is that the condensates which are generated in the RG 
flow can be treated in a simple manner. In fact, the con- 
densates appear as non-trivial saddle points in the bare 
functional integral which can easily be detected and han- 
dled by expanding around the maximum of the integrand 
[33 |. On the contrary, one always finds a convex effec- 
tive action if the effects of the possible non-trivial saddle 
points are correctly incorporated during the Legendre- 
transformation, consequently, the Maxwell construction 
hides a large part of the dynamics generated by them. 
However, the RG methods based on the effective average 
action which interpolates between the bare action and 
the full quantum effective action provide generally an RG 
flow exhibiting singularity in a truncated functional sub- 
space. The infrared (IR) singularity of the functional RG 
equation is supposed to be related to the convexity of the 
effective action for theories within a phase of spontaneous 
symmetry breaking (3^ . 



1. Wegner-Houghton RG 

The blocking in momentum space, i.e. the integration 
over the field fluctuations with momenta of the magni- 
tude between the UV scale A and zero is performed in 
successive blocking steps over infinitesimal momentum 
intervals k ^ k — Ak each of which consists of the split- 
ting the field variable, (/) — ip + (p' in such a manner that 
ip and (j)' contain Fourier modes with \p\ < k — Ak and 
k — Ak < \p\ < fc, respectively and the integration over (/)' 



leads to the Wegner-Houghton (WH) RG equation [2S 



(2 + kdk) Vk{cl)) = In (l + V^'{cl>) 



(Al) 



with Vj!'(0) = d'^Vk{4>) for the dimcnsionless local poten- 
tial Vk — k^^Vk tor d — 2 dimensions in the leading order 
of the derivative expansion, in the LPA when (j) reduces to 
a constant. (Below we suppress the notation of the field- 
dependence of the local potential.) The differentiation 
with respect to the field variable and the multiplication 
with 1 + V^' leads to the derivative form of the WH-RG 
equation, 



{2+kdk)n^-vi'{2+kdk)n- 



At: " 



(A2) 



This equation is obtained by assuming the absence of in- 
stabilities for the modes around the gliding cutoff k. The 
WH-RG scheme which uses the sharp gliding cutoff k can 
also account for the spinodal instability, which appears 
when the restoring force acting on the field fluctuations 
to be eliminated vanishes, 1 -I- Vl'{(j)) = at some finite 
scale fcsi and the resulting condensate generates tree-level 
contributions to the evolution equation. The saddle point 
(/)[) for the single blocking step k k — Ak is obtained by 
minimizing the action, Sk-Ak[(t>] = rnm^/^ (Sk[(t> + (Pq])- 
The restriction of the space of saddle-point configura- 
tions to that of the plane waves (pQ = pcos{kix) gives 

M 



Vfc-Afe(0) 



mm 

p 



duVk{(l) + 2pcos{TTu)) 



(A3) 



in LPA, where the minimum is sought for the amplitude 
p only. It was shown that the tree-level RG equation 
^3p leads to the local potential 



Vk 



1 



which can also be obtained as the solution of 
1 -I- V^'iP) = 0. 



(A4) 



(A5) 



If SI occurs during the RG flow at some scale kgi > 0, 
then Eqs. (jAl|l and (jA2|l can be applied only for scales 
k > fcgi, and the tree- level renormalization should be 
performed at scales k < fcgi. The right hand side of Eq. 
(|Aip develops a singularity when the SI occurs, but Eq. 
(|A2p does not develop such a singularity and its solution 
mathematically extends to fc — > 0. It is interesting to 
notice, that Eq. (jA2p yields the fixed-point equation 



2v; -f [v:f + j^v: ^ ci 



(A6) 



with the arbitrary constant ci, exhibiting the trivial so- 
lution V^ = ci/2 (Gaussian fixed point) and 



V, = 



1 



1 



1 



2" Svr 2 
which is equivalent to the fixed point potential ( IA4I 



(A7) 
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2. Polchinski's RG 

In Polchinski's RG (P-RG) method [29] the realiza- 
tion of the differential RG transformations is based on a 
non-linear generalization of the blocking procedure us- 
ing a smooth momentum cutoff. In the infinitesimal 
blocking step the field variable <j) is split again into the 
sum of a slowly oscillating IR and a fast oscillating UV 
components, but both fields contain now low- and high- 
frequency modes, as well, due to the smoothness of the 
cutoff. Above the moving momentum scale k the propa- 
gator for the IR component is suppressed by a properly 
chosen smooth regulator function K{y) with y = /k^ , 
K{y) -> if ?/ >> 1, and K{y) ^ 1 ff y << 1. The 
P-RG equation in LPA for d = 2 dimensions reads as 

(2 + kdk)Vk = -[n?K^ + n'h, (A8) 

where K' = dyK{y), K'^ = dyK{y)\y=o and h = 
(l/47r) K'{y). The parameters K'^ and I2 can be 
eliminated by the rescaling of the potential and the 
field variable, consequently, they do not influence the 
physics. In order to make the comparison of the RG 
flows obtained by various RG methods straightforward, 
we choose /t = — x- and Kn = —1 for which the lin- 
earized forms of Eq. (IA1[) and Eq. (jASp and the UV scal- 
ing laws obtained by WH-RG and P-RG are identical. 
Then the differentiation of both sides of Eq. (jASP with 
respect to the field variable </> yields 

(2 H- kdk)Vl = 2V;:n - ^ V;!' (A9) 

being independent of the regulator function K{y) and 
differing of the WH-RG equation (jA2p by the term 
—Vl!kdkVl with opposite sign for the non-hnear term. 

Let us note on the one hand, that the P-RG method 
treats all quantum fluctuations below and above the scale 
k on the same footing. Therefore, even if there occurs 
a scale fcsi at which 1 -I- V^" exhibits zeros, one cannot 
decide unambiguously at what scale should one turn to 
tree-level renormalization. On the other hand, Eq. (|A8[) 
with the choice of the parameters ifp — ^2 — — 
leads to the fixed-point equation similar to that of (|A6|) 
with ci = 

2v; - [vlY + -^K" = (Aio) 

47r 

and exhibits trivial fixed-point solutions: the Gaussian 
one (14 = const) and the high-temperature (or infinitely 
massive) fixed point 

K = - ^ (AH) 

which is only accounted for the P-RG method. Let us 
note, that (|A11[) is similar to (|A7[) which is obtained in 
the framework of the WH-RG method but with opposite 
sign. For dimensions d = 2 it was shown in [s^ that non- 
trivial fixed point solutions in LPA are either singular at 
finite if or periodic. 



3. Effective average action RG 

The effective average action (EAA) RG method [Sll, 
m, [3^ has grown out of the idea of coarse-graining the 
quantum fields and it interpolates between the bare ac- 
tion and the full quantum effective action. The scale- 
dependent effective average action Ffc satisfies the func- 
tional differential equation 

kdkVk = ^Tr (r*^'^ + i?fe) kdkRk (A12) 
(2) 

where denotes the second functional derivative of 
the effective action. Here Rk is a properly chosen IR 
regulator function which fulfills a few basic constraints 
to ensure that F^ approaches the bare action in the UV 
limit (fc — > A) and the full quantum effective action in the 
IR limit (fc ^ 0) and to guarantee that no IR divergences 
are encountered in the presence of massless modes. In 
the present paper we shall use Rk{p^) = p^rOy) with 
the power-law regulator r{y) = y~'' {b > 1) [33| and 

the optimized regulator [1] r{y) — — 1^ 6(1 — y) (9 

denotes the Heaviside step- function) where y — p^/k"^. 
For d — 2 dimensions Eq. (jA12[l can be rewritten in the 
LPA as 

{2 + kdk)V,^~— dy— (A13) 

47r Jo y(l + r) + V^' 

with r' = dr/dy for the dimensionless local potential. 
The truncation of the basis set of functions in which the 
local potential is expanded may introduce a dependence 
of the physical results (e.g. the critical exponents) on the 
particular choice of the regulator function. In the litera- 
ture several optimization procedures have been proposed 
in order to achieve better convergence of the critical ex- 
ponents with the removal of that truncation e.g., Refs. 

For example, if one compares various RG schemes, 
it has been argued that the fastest convergence can be 
achieved by using the optimized regulator [6]. Concern- 
ing the power-law regulator it was argued that the opti- 
mal choice of the parameter b is met when the minimum 
value 

^(^) = (6-l)Vi)/^ (A14) 

is maximal, which happens for 6 = 2, i.e., for the quartic 
regulator. For arbitrary parameter value &, the propaga- 
tor in the right hand side of Eq. (jAlSp may develop a 
pole at some scale ksi and at some value of the field <j) for 
which Vl'{(j>) = —C(b) holds, which signals the occurring 
of SI. It was shown that in such a case one has to seek the 
local potential for k < ksi by minimizing Ffc in the sub- 
space of inhomogeneous (soliton like) field configurations 
and ends up with the result [33 ] 

Vk = -^C(6)02 (A15) 
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of parabolic shape, the solution of the equation 

Cib) + 14" (0) - 0. (A16) 

For the quartic regulator the integral over y can be 
performed and Eq. (jA13p rewritten in the explicit form 



(2 + kdk)Vk 



arctg 



iln 



1 



2W|l-(^n")'l 



(A17) 



for 



for (iy," 



< 1 



> 1 



when the limit A — > 00 is taken. Making use of 
I'Hospital's rule it is straightforward to show that the 
right hand side of (|A17|) tends to for iy^" —> 1 from 
either below or above and, consequently, the couplings 
behave non-analytically at such scale k^, but do not di- 
verge to infinity. However, if there exists a scale fc_ > 
for which \VIJ ^ -I, then the right hand side of (|A17|) 
becomes infinite which indicates the appearance of SI 
due to vanishing of the inverse of the IR regulated prop- 
agator. Therefore the scale fc_ can be identified with 
the scale of the SI, fc_ = k^i- The scale fc+ may exist 
even if no SI arises during the flow, it is an artifact of 
the regulator used, but of no physical significance. We 
shall demonstrate this latter case for the one-component 
scalar field theory with polynomial interaction. The two- 
dimensional generalized sine-Gordon model shall demon- 
strate another interesting case when the scales k± are 
equal. 

It is worthwhile noticing that Eq. (|A13p with the 
power-law regulator leads to the WH-RG equation (jAl[) . 
and Eq. (jA16P leads to Eq. ()A5j) for 6 = 1 as well as for 
6 — !■ 00 in the limit A ^ 00. This feature of the case 
h ~ \ holds only for d = 2. 

The RG equation with the optimized regulator in the 
LPA reads as 



(2 + kdk)Vk 



1 



47r 1 



(AI8) 



where a field-independent constant has been added to 
the right hand side of the equation. Let us note that 
if SI occurs, the right hand side of Eq. (jAlSp becomes 
infinite at the scale fcsi, so that all or some of the cou- 
plings should tend to infinity at this scale. One should 
again turn to the explicit treatment of SI. The equation 
obtained by deriving (jAlSp with respect to the field and 
multiplying its both sides by {y + Vl!)'^ leads to the fixed- 
point equation 2(1 Vl'fVi = which exhibits 
the trivial solution \4 — const (Gaussian fixed point), 
and also the parabolic solution (|A7p . Let us note, it was 
demonstrated that Eqs. (jA18|) and (jASp give the same 
critical exponents for the 0(N) symmetric scalar theory 
in three dimensions and the two RG equations can 



be transformed onto each other via a suitable Legendre 
transformation However, the RG trajectories and the 
singularity structures of (jAlSP and (|A8P could be differ- 
ent. 



4. Functional Callan-Symanzik RG 

In the functional Callan-Symanzik (CS-RG) type in- 
ternal space RG method [soj, the successive elimination 
of the field fluctuations is performed in the space of the 
field variable (internal space) as opposed to the usual RG 
methods where the blocking transformations are realized 
in either the momentum or the real (external) space. This 
can be achieved by introducing an additional mass term 
into the bare action, 



(A19) 



with the control parameter A. For A = Aq being of the 
order of the UV cutoff the large-amplitude fluctuations 
are suppressed and decreasing the control parameter A 
towards zero, they are gradually accounted for. The func- 
tional evolution equation for the effective action is 



1. 



-Tr 



(2) 



2A^ 



(A20) 



where F^^ ^ = x/ SipS^p. Eq. (|A20p is equivalent to 
the RG equation (jA12p for the EAA with the power- 
law regulator function with 6=1 for X — k. By using 
Eq. (|A13p the functional Callan-Symanzik RG equation 
for the one-component scalar fleld theory for dimensions 
d = 2 in the LPA reads 



{2 + Xdx) Vx = --Lln 



1 + Vi' 



(A21) 



where V\ = X ^Va is the dimensionless scale-dependent 
effective potential. Adding the fleld-independent term 
-l/(47r)ln(A2/A2) to the right hand side of Eq. (|MT|l . 
one can take the limit A ^ 00, 



(2-f A^a) Vx 



■In 1 



(A22) 



This equation is mathematically equivalent to the two- 
dimensional WH-RG equation in the LPA assuming the 
equivalence of the scales A = fc. However, for dimensions 
d ^ 2 the functional Callan-Symanzik RG and the WH- 
RG differ from each other. Assuming the above men- 
tioned equivalence of the scales A and fc, there occurs the 
same singularity in the right hand side of (jA22p as the 
one in the WH-RG approach. Therefore, the functional 
Callan-Symanzik RG signals the SI with the vanishing of 
the argument of the logarithm in the right hand side of 
(|A22p . The solution of (|A22p provides the scaHng laws 
down to the scale fcsi and one has to turn to the tree-level 
renormalization with the help of the WH-RG approach 
in order to determine the IR scaling laws. 
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APPENDIX B: COMPARISON OF RG SCHEMES 
FOR POLYNOMIAL INTERACTION 

Here, we compare the applicability of the frequently 
used functional RG methods to two-dimensional one- 
component Euclidean scalar fields with polynomial self- 
interaction. We show that the flow of the polynomial 
model determined by the P-RG differs of that obtained 
by other RG methods in the deep IR regime. For the 
polynomial local potential 

^ 1 

^'=('^)=E 7^32n(fc)0'" (Bl) 
n=l ^ 

the scale-dependence is encoded in the dimensionless 
coupling constants g2n{k), related to their dimensionful 
counterparts via k^g2n — g2n- The bare action exhibits 
the Z2 symmetry — > — (/>, and it is well-known that there 
exist two phases of such a model with either unbroken or 
broken Z2 symmetry. 

Let us start with the discussion of the RG flow in the 
phase with unbroken Z2 symmetry. Inserting the po- 
tential dHU into any of the RG equations (fA2|l . (|A18|1 . 
(|A17[) . (|A9[) and expanding both sides of them into 
Taylor-series in the dimensionless field-variable 0, one 
arrives at a coupled set of ordinary, non-linear differ- 
ential equations for the couplings g2n{k). These have 
been solved numerically for the same initial condition 
52(A) = -0.001, 54(A) ^ 0.01, g2n>A = at the UV cut- 
off A which ensured that no SI occurred, i.e. the inequal- 
ity + g2{k) > has been kept during the flow. The 
RG flow of the couplings 32 (fc) and g4,{k) for various RG 
methods and for iV = 10 is plotted in Fig.[8l There were 
no appreciable changes in the results when we increased 
N further. The WH-RG, the EAA-RG with optimized 
regulator, and the EAA-RG with quartic regulator give 
qualitatively the same results, the dimensionful couplings 
become constant in the IR limit. Let us note that g4^{k) 
is non-analytic at the scale for which k^ — g2{k+) if the 
flow is determined by the quartic regulator RG. 

For the P-RG method the flow of the couplings in the 
IR region (fc 0) differs of that obtained by WH-RG 
and EAA-RG methods (see Fig.jS]). The RG trajectories 
tend to the trivial high-temperature (or infinitely mas- 
sive) fixed point (|Alip which is accounted only for the P- 
RG method. (Let us note, that in three-dimensions the 
critical behavior is dominated by the non-trivial Wilson- 
Fisher fixed point.) However, one may argue that the RG 
analysis is physically uninteresting in the IR limit since 
below the momentum scale fc^ = 52 (fc) the mass term 
suppresses the quantum fiuctuations in the propagator 
and after a rather small transient domain, the RG fiow 
obtained by the other RG methods becomes trivial, i.e. 
all the dimensionless couplings (52, 54) scale as ~ fc^^ (for 
d = 2), consequently, the corresponding dimensionful pa- 
rameters (52, 94) tend to constants, see Fig. [8] It is also 
demonstrated in Fig. [51 that even the P-RG shows up the 
change of the sign of 52 (fc) because it occurs above the 



mass-scale. However, in this paper we show that there 
may be special situations such as for example the case of 
the MSG model where one has to go beyond the mass- 
scale in order to map out the phase structure of the model 
in a reliable manner. In such a case the P-RG method 
is inappropriate for quantitative analysis. 




FIG. 8: The RG flow of the dimensionful couplings g2{k) and 
(74 (fc) obtained by various RG methods for the same initial 
condition g2{A) = -0.001, sf4(A) = 0.01, S(2n>4(A) = 0. The 
full, dashed, and dashed-dotted lines correspond to the EAA- 
RG with quartic regulator, the EAA-RG with optimized reg- 
ulator, and the WH-RG flow, respectively and the dotted 
lines represent the flow obtained by the P-RG method. The 
mass g2{k) changes its sign in any cases at (fc/A) ~ 3T0~^ for 
the given initial conditions. If the flow is determined by the 
quartic regulator RG, Qiik) is non-analytic at the momentum 
scale k+ ~ 2.64 ■ lO"'^ for which fc^ = g2{k+). 

The discussion above has been restricted to the phase 
with unbroken Z2 symmetry, when no SI occurs during 
the RG flow. The RG flow has also been investigated in 
the literature for the phase with spontaneously broken 
Z2 symmetry in the frameworks of both the WH-RG 
involving tree-level renormalization and the EAA-RG. 
The truncated WH-RG flow becomes numerically unsta- 
ble, i.e. the couphngs g2n{k) start to heavily oscillate 
approaching the scale fcsi from above where the SI oc- 
curs. But it has been shown that the tree-level evolution 
for RG trajectories started below the scale fcgi run into 
the universal IR fixed-point potential (IA4|) . Applying the 
explicit treatment of SI proposed in 32], one would ob- 
tain similar results by the EAA-RG. It has been argued 
in 32] that the dimensionless IR effective potential for 
k < kgi reads as the fixed-point potential (jA15l) (c.f. Ap- 
pendix lA 3|) . It was also demonstrated in [38] that the 
appearance of SI can be avoided by a suitable rescaling 
of the RG equations at least in case of the polynomial 
scalar theory. However, an attractive IR fixed point ap- 
pears in the rescaled RG flow which can be identified as 
the effective potential (jA15l) . 
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